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Transport measurements allow sensitive detection of nanomechanical motion of suspended carbon 
nanotubes. It has been predicted that when the electro-mechanical coupling is sufficiently large a 
bistability with a current blockade appears. Unambiguous observation of this transition by current 
measurements may be difficult. Instead, we investigate the mechanical response of the system, 
namely the displacement spectral function; the linear response to a driving; and the ring-down 
behavior. We ffnd that by increasing the electro-mechanical coupling the peak in the spectral 
function broadens and shifts at low frequencies while the oscillator dephasing time shortens. These 
effects are maximum at the transition where non-linearities dominate the dynamics. These strong 
signatures open the way to detect the blockade transition in devices currently studied by several 
groups. 

PACS numbers: 73.23.Hk, 73.63-b, 85.85.+j 


Recently enormous progress has been achieved in the 
detection of carbon nanotubes (CNT) bending modes 
by electronic transport measurements [iHia- Since nan¬ 
otube oscillators have remarkable mechanical properties, 
devices with record mass [7] and force m sensitivity have 
been realized. Transport experiments allow information 
to be obtained on the mechanical mode by measuring 
different quantities. The main ones are the oscillation 
amplitude in response to an external drive [T], the os¬ 
cillator displacement spectral density, Sxx{^) cni, and, 
more recently, the ring-down time of the oscillator m- 

The recent experimental advances allow one to view 
the behavior of such systems in the strong coupling limit 
from a new perspective. Defining Fq as the difference of 
electrostatic force acting on the nanotube when one elec¬ 
tron is added to the suspended part and k as the spring 
constant for its displacement, one can introduce a po- 
laronic energy scale sp = F^/k. For a classical system 
resonating at pulsation ujq (much smaller than the bias 
voltage eV or the temperature T) it has been predicted 
[HHUI that, if ^ eU and T, the current in the device 
can be blocked and a bistability can appear (e is the elec¬ 
tron charge and we set both the Planck and Boltzmann 
constant to 1). The energy sp can be estimated for actual 
experiments: For a CNT of 1 /am length, 1 nm radius, 
and suspended at a distance d = 500 nm from a gate one 
finds uJol27T « 50 MHz, Fq = 10“^^ N, /c = 4 • 10“^ N/m, 
and thus sp ^ 16 mK. One should thus work at very 
low temperatures ^ 16 mK: This is probably why the 
current blockade transition has not yet been observed for 
mechanical bending modes. While a (Franck-Condon) 
blockade [ig ng has been observed m for breathing 
modes in the regime of incoherent transport. To increase 
ep one can reduce the distance of the CNT from the gate 
electrode or operate the system close to the Euler buck¬ 
ling instability ED Eg. The energy ep scales quadrat- 
ically with d, it is thus realistic to increase this energy 
up to the Kelvin range by reducing the distance <7 to 100 
nm. In any case, a clear observation of the transition 


will require temperatures of the order 100 mK. At such 
low temperatures the typical tunnelling rate F becomes 
larger than T, leading to coherent transport through the 
CNT. In this regime the current blockade can take place 
only if ep is larger than a critical value 5c of the order 
of F [T5] . Since one of the main experimental diffi¬ 
culties is to reach large values of 5p, the case 5p ^ F is 
particularly interesting. The transition could then be in¬ 
vestigated at fixed and low V and T by varying ep, that 
can be tuned with the gate voltage. One drawback of 
this limit is the large width F of the electronic level: The 
conductance dependence on 5p is smooth on the scale F 
failing to provide a clear indication of the transition. 

In this Letter we show that, similar to critical phenom¬ 
ena, the transition can be better investigated by looking 
at the behavior of the phonon mode that becomes soft for 
5p = 5c. We study Sxx{^)^ Itie driving response function, 
and the dephasing and ring-down time as a function of 
the coupling constant ep. We find that all of these quan¬ 
tities have a very peculiar behavior at the transition. The 
dynamics of the mechanical mode is dominated by non¬ 
linear terms leading to a separation of time scales that 
is maximal at the transition. The theory presented gives 
clear indication on how to unambiguously observe the 
transition using available methods of measurement. 

The model. We consider a suspended CNT. We assume 
that a single electronic level is relevant for transport. We 
neglect the spin degrees of freedom. The Hamiltonian 
reads: 

2 7 

H = HL+HR + HT + {eo-xFo)Sd + + W > W 

where d is the destruction operator for the electronic 
level on the dot, x is the displacement of the relevant 
mechanical mode, p the conjugated momentum, m the 
mode effective mass, k the spring constant (giving a 
pulsation cjq = \/klm) and Fq the electrostatic force 
acting on the dot when an electron is added. The 
first three terms describe the leads and their coupling: 
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Ha = Y^ki^ak - with a = L and R, for left 

and right lead, Cak the electronic spectrum, fio, the chem¬ 
ical potential; and Ht = +h.c. the tunnelling 

Hamiltonian. From these quantities one can define the 
single-level width with the density of 

states and F = Tl -\-Tji. 

In the Born-Oppenheimer limit (Ta ^ c^o) the dis¬ 
placement of the mechanical mode can be described by a 
Langevin equation: 

mx + A{x)x + mujQX = Fe{x) -h ^{t) , (2) 

where the dissipation A{x)^ the average force Fe{x) = 
and the stochastic force ^(t) are due to the 
electrons tunneling through the quantum dot Hi US]. 
The explicit expressions for H, Fe, and = 

D{x)6{t — t') have been obtained in Ref. [15]: 
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( 3 ) 


A{x) = - eo - Foxf +T^]^, and 

D{x) = FSrLTii/iTTr^)[h{nL) - Hur)], where h{n) = 
arctan 2 : -h z/{z‘^ + 1) with z = {p — eo — Fox)/r. In the 
same limit a Fokker-Planck equation for the probability 
P(x,p, t) can be derived [23l|^: 

dtP = - FdpP + -dp{pP) + ^dlP , (4) 

m m 2 ^ 

with F{x) = Fe{x) — kx. 

Softening of the mechanical mode. We assume that 
the device is symmetric: Fq, = F/2. The presence of the 
mechanical coupling modifies the electron-hole symme¬ 
try point for eo to the value eo = {pl F Pr)I‘^ F SpI2. 
We will always assume this value from this point on. 
Defining y = x — Fo/2k^ one determines F{y) = —kyF 
(Fo/27r) arctan[(Fo^ + aeV/2)/T] which depends 

only on the bias voltage eV = Pl — Pr and is anti¬ 
symmetric in y. The equilibrium positions are defined 
by the solutions of the equation F{y) = 0. The line 
sp = s^(V) = 7rr[l + (eR/2F)2] [for eV/T < 2/^/3] sepa¬ 
rates the monostable region from the bistable region (see 
inset of Fig. 0). 

Let us now define ujrn at a stable point yg as = 

— (dFjdy^y^. It goes smoothly from ujq to 0 for 5p < 5c 
with the analytic form = [Sc — ^p) j while for 

^ it reads 2{ep-ef)lec [see Fig. 0 orange dashed 
lines]. 

The vanishing of ujm suggests that its direct measure¬ 
ment should allow the detection of the transition with 
great accuracy. As in phase transitions, this mode be¬ 
comes soft, leading to a strong response at the transition. 
The dip in the gate voltage dependence of uJm observed 
by four different groups mmani] is the precursor of this 
softening. Nevertheless, one should be cautious since the 
definition of uJm only takes into account the first deriva¬ 
tive of the force at the minimum of the potential. When 


FIG. 1. Density plot of Sxx as a function of u and sp. The 
values of cUm (orange dashed line), 2(jUm (blue dot-dashed line), 
and 3ujm (red dotted line) are shown. The units of Sxx are 
^zpm/^o = (mc^o)”^, where Xzpm = {muJo)~^^‘^ is the zero- 
point motion displacement. The symmetry of the potential 
implies that only odd harmonics are present for sp < Sc- 
Inset: phase diagram in the plane eV-Sp for the stability of 
the effective potential. 


this term vanishes the next order terms in y become im¬ 
portant and the response of the system can no longer be 
predicted simply by the value of uJm- Therefore, in the 
following we calculate the typical measurable quantities 
and study their behavior when Sp is swept through the 
transition. 

Fluctuation spectrum. We define the displacement 
fluctuation spectrum Sxx{^) = f e*^^(it(x(t)x(0)), with 
x{t) = x{t) — {x). This quantity has been measured re¬ 
cently in Ref. m- We can obtain Sxx numerically from 
the Fokker-Planck description following the method used 
in Ref. [17]. Writing Eq. ^ 8iS dfP = jCqP the spectrum 
takes the form: 


F,,(c^) = -2Tr 


Co 


d^Fst 


( 5 ) 


where P^t is the stationary solution of the problem, sat¬ 
isfying both Co Pst = 0 and the normalization condition 
TrPst = 1. The operator x is defined as xP = xP{x^p). 

Let’s begin by discussing the stationary solution of 
Eq. Psf In agreement with similar models [22l 

|25] we find that for sufficiently small eR, even if the 
system is out of equilibrium, the stationary distribu¬ 
tion function takes the simple Gibbs form Pst(x,p) = 
A/'expl—P(x,p)/Teff}, with Teff = eR/d, JV a normaliza¬ 
tion factor, P{x,p) = p^12m + t/(x), U{x) = —dP/dx^ 
and U{x) = 0 at its minimum. This result is due to the 
smooth dependence on x of both D{x) and A{x) on the 
scale of the spread of the probability distribution P(x,p) 
for eV F. 

We come now to the displacement spectrum obtained 
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from Eq. ^ that we show for eVjV — 5 • 10“^ in Fig. 
and 1^ As anticipated, as coupling increases the reso¬ 
nance broadens and shifts at low frequency. More sur¬ 
prisingly, the peak position shows a minimum at a small 
but finite value of cj at the transition with a maximal 
broadening (cf. Fig. |^b). The position of the minimum 
and its width depend on the bias voltage. A strong tele¬ 
graph noise appears for 5p > 5c (dark region at cj ^ 0), 
signaling the hopping of the systems between the two 
minima in the potential [I3II7I. This is also a strong in¬ 
dication of the transition [see Fig. Q-d]. In the bistable 
phase [sp > Sc) a double peak in the spectral function is 
visible. 

In order to understand this behavior we take advan¬ 
tage of the separation of time scales of the problem. The 
damping A{x) and the fluctuations D{x) are both gen¬ 
erated by the non-equilibrium electronic transport and, 
by hypothesis, are parametrically smaller (as uJo/V) than 
the Hamiltonian terms in Eq. This implies that the 
system performs many oscillations on the closed trajec¬ 
tory in the phase space that satisfies E{x{t),p{t)) = F, 
before drifting to a nearby trajectory on the slow time 
scale where yp = f dxP^t{x^p)A{x)/m is the av¬ 
erage dissipation coefficient [mils]. For each energy E 
one can then calculate the pulsation of the closed trajec¬ 
tory u;(F) = 1/[27r(m/2)^/^ . The 

non-linearities present in U(x) induce dispersion in cj(E). 
Then, from our definition, cj^n indicates only co’(O). Ener¬ 
gies up to eV are populated as an effect of the stochastic 
fluctuations. They all contribute to the fluctuation spec¬ 
trum leading to an inhomogeneous broadening that spans 
the frequencies between cj^n and uj(eV). 

In order to provide a quantitative verification of 
this interpretation we calculate the spectrum by ne¬ 
glecting the effect of the dissipation and consider¬ 
ing only the interference of the different trajecto¬ 
ries populated according to Pgt [27]. This gives 
Sxx(t) = f dxodyo-Pst(xo,Po)x(t)x(0), where x(t) satis- 
fies the equation of motion mx = F(x), with the initial 
conditions x(0) = xq, and j:(0) = p^jvn. In terms of 
the Fourier coefficients of fixed energy periodic trajecto¬ 
ries [xpit) = the spectrum takes the 

form: 

pOO 

SU^)= P{E)dEy22TTS{uj-nuj{E))xl{E), (6) 

with P{E) = Afe-^/^^<‘2TT/!x{E). 

Two limits can be analyzed [58]: When the quar- 
tic term is much smaller than the quadratic one, the 
full width at half height of the resonance is Aoo ^ 
l.buJo{eV/Sc){sp/Sc)‘^{l — 5p/5c)“^/^ with a small posi¬ 
tive shift of the maximum from of the same order. In 
the opposite limit of vanishing harmonic term {sp = Sc) 
the potential can be approximated as quartic. In this 
case Acj/cjo = 0.50(eF/r)^/^ with a peak position ujm 
at 0.85c(;o(eF/r)^/^. This voltage dependence can be re¬ 
lated to the dispersion of uj{E) that vanishes as E^^^ for 



FIG. 2. Comparison of the full numerical solution of the 
Fokker-Planck equation for Sxx (blue solid lines) with the 
one obtained with Eq. § (dots) for spIsc = 0.87, 1, 1.11 
(a, b, and c panel, respectively) in units of (mc^o)”^. The 
numbers label the order of the harmonic, while the letters h 
and 1 in the c panel indicate the high- and low-frequency con¬ 
tributions. Panel d: *S'cccc(0) as a function of ep/F indicating 
the onset of the telegraph noise at the transition. 


ep = Sc- Remarkably, at criticality the Q-factor of the os¬ 
cillator defined as ujm/A uj takes the universal value 1.71, 
independently of F or F. The crossover between the two 
regimes takes place for 1 — 5p/5c ~ 1.71(eF/5c)^/^, thus 
for the values considered in Fig. [^the quartic region is 
restricted to l — spjec < 0.03. Finally, the double peak of 
the spectral function for 5p > 5c can be explained by the 
contributions of the low-energy high-frequency trajecto¬ 
ries around each single minimum, and those, at higher 
energy and lower frequency, revolving around both min¬ 
ima. The comparison with the numerical result presented 
in Fig. (1) shows a very good agreement. 

Driving. Let us consider the other main tool used 
to detect mechanical motion: The response to a driving 
force of frequency ujp. We can find the linear response 
of the system by letting E{x) E{x) + Fp cos(co’pt) 
in Eq. The evolution operator becomes C{t) = 

Cq -\- 2Cd cos(ci;pt), with Cp = —Epdp/2. After a tran¬ 
sient time the solution can be written as a Fourier se¬ 
ries P{t) = e'^^^^^Pn where each Fourier component 

can in turn be expanded as a power series of the driv¬ 
ing parameter Ep: P^ = YX={)Pn,ki with Pn^k of or¬ 
der Fh- This leads to the equation for each component 
{inujp — £o)^n,o = 0 and 


{inujp - £o)^n,/c+l = ^D{Pn+l,k + (7) 


with the condition TrP(t) = 1 . Eq. 0 can be solved by 
recursion. The time dependence of the displacement then 


reads x{t) = Tr 


xP{t) 


Edx{^d)P^^^ + C.C., where 


X(w)=Tr xPiit) 


= Tr 


x(^icj Eq') ^pFst 


(8) 
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FIG. 3. Time dependence of {x{t)) (green solid line) and its 
envelope (blue solid line) for epjsc = 0.32, 1, and 1.11 (a, b, 
and c panel, respectively). The exponential decay on the scale 
2/'yE is shown red dashed. In panel a and b the orange dotted 
line gives the result of the analytical expressions discussed in 
the text. Panel d: comparison of the £p-dependence of Acj, 
27 c^, and yp. 


Naturally, the relation between x(co’) and Sxx{^) comes 
into question. If Pst has a Gibbs form then FjjCqxP^i = 
—2Teff>Ci:)Pst. This leads to a fluctuation-dissipation re¬ 
lation: 


Im [x(w)] = . (9) 

Thus, for eV <C F, x and Sxx give access to the same 
information in two independent ways. For larger voltages 
expression ^ always holds while Eq. ^ will be violated. 

Ring-down behavior. Finally, let us consider the re¬ 
sponse for time t > 0 of the oscillator when the coherent 
drive is switched off at t = 0. A damped harmonic oscil¬ 
lator relaxes exponentially on a time scale (the ring-down 
time) given by the same dissipation coefficient that also 
determines the width of the resonance of the response 
function. For nano-mechanical oscillators it has recently 
been shown m that this may not be the case. Non- 
linearities induce frequency noise, which in turn is re¬ 
sponsible for phase fluctuations of x{t). The average over 
many realizations of x{t) decays then on the time scale 


7 “^, the value for which phase fluctuations become the 
order of 27r. Since the energy is insensitive to the phase, 
its average decays on the same time scale as the single 
realization 7 ^^. 

With respect to our problem, we use the solution of 
the Fokker-Planck equation with driving [Eq. 01 as the 
initial condition and then find x{t) and E{t) from the 
evolution of the probability with jCq. The result as a 
Laplace transform reads: 


{x{z)) = Tr x{z — jCq) ^P{t = 0) . (10) 

In a similar way we can calculate also the Laplace trans¬ 
form of the evolution of the total energy by letting x 
E{x^p) in Eq. ( p^ . One can then obtain the time depen¬ 
dence by numerically implementing the Cauchy theorem 
{x{t)) = f^{x{z))e~^'‘^dzl{27ri), where (7 is a contour that 
encloses the poles of {x{z)) for Re^ < 0 . 

We find that the energy exponentially decays on the 
scale 7 ^^, even at the transition. On the other hand, as 
shown in Eig. 0, {x{t)) decays on a much shorter scale 


that we define 7 ^^- Eig. (|^-d shows the ep dependence 
of Acj, 27 ^, and 7 ^;. The width Acj, obtained from the 
form of Sxx^ coincides within the numerical accuracy with 
27 ^, proving that frequency noise is the responsible of 
the faster decay of {x{t)). Both present a pronounced 
maximum at Sp = 5c, indicating the transition. Using 
the approach presented for the analytical calculation of 
Sxx{^) we find that {x{t)) decays as 1 /( 1 + 7 ^)^, where 

7 c^ = O.dlAcj [29]. Similarly for 5 p = 5c the decay scale 
is proportional to Auj with the analytical form given by 
the dotted line of Eig. 

Conclusions. We found that the study of the me¬ 
chanical properties of the suspended carbon nanotube 
open new perspectives for the observation of the current- 
blockade transition occurring at low temperaure and volt¬ 
age for 5 p = ttF. Indeed, for that value the quadratic 
part of the effective potential vanishes, leading to strong 
frequency and phase fluctuations that remarkably mod¬ 
ify the typically measured response functions {Sxx^ X? 
and ring-down behavior). The Q-factor of the resonator 
takes the minimum and universal value 1.71 at critical¬ 
ity, where the separation of time scales is also maxi¬ 
mal. These results can lead to the observation of the 
current-blockade transition by mechanical measurements 
in devices currently investigated by several experimental 
groups. 
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